Effect of plasma resonances on dynamic characteristics of double graphene- layer 

optical modulator 
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We analyze the dynamic operation of an optical modulator based on double graphene-layer(GL) 
structure utilizing the variation of the GL absorption due to the electrically controlled Pauli block- 
ing effect. The developed device model yields the dependences of the modulation depth on the 
control voltage and the modulation frequency. The excitation of plasma oscillations in double- 
GL structure can result in the resonant increase of the modulation depth, when the modulation 
frequency approaches the plasma frequency, which corresponds to the terahertz frequency for the 
typical parameter values. 



I. INTRODUCTION 



The gapless energy spectrum of graphene layers 
(GLs)) [if, results in the interband absorption of the elec- 
tromagnetic radiation from the terahertz to ultraviolet 
range [2j. This opens up prospects to use graphene struc- 
tures in active and passive optoelectronic devices. Novel 
lasers, photodetectors, modulators, and mixers have been 
proposed and studied [HIH (see also the review pa- 
per (l6j and references therein). Varying the controlling 
voltage one can effectively change the electron and hole 
densities and the Fermi energy and, hence, the intraband 
and interband absorption in GLs. An increase in the elec- 
tron (hole) density increases the intraband absorption 
(associated with the Drude mechanism), but it decreases 
the interband absorption due to the Pauli blocking ef- 
fect. Which mechanism dominates depends on the in- 
cident photon energy Ml and the momentum relaxation 
time of electrons and holes, r. At the frequencies of in- 
frared and visible radiation, f2 3> r _1 , so that the Drude 
absorption is weak. Even in the terahertz (THz) range, 
the latter inequality can be valid for high quality GLs 
(like those studied in Refs. [13, 13)) even at room tem- 
peratures. 



'Electronic mail: v-ryzhii(at)riec. tohoku. ac.jp 



In this paper, we develop a device model of a double- 
GL optical modulator proposed and demonstrated in 
Ref. . In this modulator, GLs were separated by 
relatively thick barrier and the structure was integrated 
with an optical waveguide The operation of the device 
under consideration is associated with the filling of GLs 
with electrons and holes injected from the contacts un- 
der the self-consistent electric field created by the ap- 
plied voltage and the electron and hole charges in GLs. 
This process determines both static and dynamic char- 
acteristics of the double-GL modulator. At the non- 
stationary conditions, the dynamics of the electron-hole 
plasma in double-GL can exhibit resonant response due 
to the excitation of plasma oscillations similar to those 
well known in more traditional two-dimensional electron 
and hole systems |20l43C| . The plasma oscillations in GL- 
structures also were considered previously (see, for in- 
stance, Refs. (3ll-l35|V A resonant plasmonic THz using 
double-GL structures was recently studied in Ref. [36j . 

In this paper, we develop a device model for double-GL 
modulators of optical radiation and demonstrate that the 
resonant excitation of plasma oscillations in double-GL 
modulators can provide an efficient modulation of optical 
radiation by high frequency signals, in particular, in the 
THz range. 

The variation of the Fermi energy in double-GL is asso- 
ciated with the electron and hole injection and extraction 
by the side contacts. Therefore, the consideration of the 
electron- hole plasma dynamics in double- GLs must ac- 
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FIG. 1: Schematic view of (a) double-GL modulator struc- 
ture coupled with an optical waveguide and (b) band diagram 
of GLs under a voltage drop between them (shaded areas in- 
dicate the states occupied by electrons) . 



count for the self-consistent electric field found from the 
solution of the hydrodynamic equations and the Poisson 
equation. 
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where c is the speed of light in vacuum and k is the dielec- 
tric constant of the waveguide material. The symbol (...) 
means the averaging accounting for the distribution of 
the optical field En(x,y) in the waveguide, where the x- 
axis and y-axis correspond to the direction along double- 
GL structure and perpendicular to the z-axis correspond- 
ing to the direction of the wave propagation direction in 
the waveguide. Thus, 



(CTfi)rn = 



f_ L Reaa\E Q (x,0)\ 2 dx 
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where 



is the mode overlap factor and 2L is the GL length, which 
is approximately equal to the spacing between the side 
contacts as shown in Fig. 1(a). 

We assume that the electron and hole momentum re- 
laxation time is associated with the scattering due to 
disorder and acoustic phonons, so that its energy depen- 
dence is given by r _1 = v(e/To), where v is the char- 
acteristic scattering frequency at V — (at the Dirac 
point). In this case, the real part of the conductivity of 
two GLs at SI z^o can be presented as 1Q| 
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II. EQUATIONS OF THE MODEL 



We consider the double-GL modulator reported in 
Ref. Its simplified structure is shown in Fig. 1(a). 
We assume that each GL is connected to one side con- 
tact and is isolated from the opposite contact (con- 
nected to the other GL), The voltage, V mi applied be- 
tween these contacts, V m = Vo + 5V m , where Vq and 
SV m are the bias and modulation components. We set 
SV m (t) = 5V m exp(— iuit) 1 where lj is the modulation fre- 
quency. The latter is much smaller than the frequency 
of the incident optical radiation £7: uj <C O. The system 
of two highly conducting side contacts can be considered 
as a slot line enabling the propagation of the modulation 
signals. The side contacts can also be connected to or 
be a part of THz antenna, which converts the incoming 
THz radiation into the modulation voltage. 

The absorption coefficient of the light wave propagat- 
ing along the waveguide with double-L on topGL is de- 
termined by the real part of the double-GL conductivity 
Recrn. 
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Here /i+ and \i- are the GL Fermi energies counted from 
the Dirac point in the upper and lower GLs, respectively, 
T is the temperature, and e is the electron charge. The 
first five terms in Eq. (4) correspond to the contributions 
of the interband transitions to the conductivity of the up- 
per and lower GLs. These terms explicitly account for the 
Pauli blocking effect. The last term in Eq. (3) accounts 
for the intraband transitions. It is presented in the form 
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providing an interpolated dependence of the intraband 
conductivity on the Fermi energies and the temperature. 
At slow variation of the applied voltage, /u+ = //_ = fj,, 
where fi obeys the following equation: 



2//- 



eVrr, 



(4) 



The quantity A is determined by the electric field be- 
tween GLs and the thickness of the barrier layer between 
GLs d. Generally speaking (at sufficiently fast variations 
of the voltage), the spatial distributions of the electron 
and hole densities do not follow the variation of the ap- 
plied voltage, and the Fermi energies depend on the co- 
ordinate x: /i + = H+(x) and = fi-(x). 

The Fermi energies in GLs are governed by the follow- 
ing equation: 



At sufficiently large bias voltage Vq, the electron and 
hole systems in GLs become degenerate (i.e., /i > T), 
and Eq. (10) yields 



2V m kV m 
so that, taking into account Eq. (4), 



A ~ < 1 "„, - 2T\/^jp = eV m - hv w \j 
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Using Eqs. (1) and (3) and considering Eq. (8), we 
arrive at the following formula for the absorption coeffi- 
cient: 
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where vw = 10 8 cm/s is the characteristic velocity of 
electrons and holes in GLs. Since the spacing, d, between 
GLs is rather small compared to 2L, one can use the 
following formulas which relate the densities E± and the 
electric potentials of GLs <p± : 
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The negative values of S± correspond to the case when 
a GL is filled by electrons, whereas the positive values 
correspond to the filling by holes. 



III. MODULATION CHARACTERISTICS 



When the control voltage V m varies slowly, one can set 



ip± = ±- 



and, hence, 



S± = T 



k Vm 
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where /x is governed by the following equation: 
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If d = 10 nm, k = 7 (A1 2 3 ), and T = 300 K, from 
Eq.(10) we obtain V ~ 30 mV. 
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Here 



2naTn 
Vk 



(13) 
(14) 



v is the characteristic collision frequency of electrons and 
holes (which is assumed to be proportional to T), and 
a = e /cH ~ 1/137 is the fine structure constant, (so 
that ira ~ 0.023). 

In the near infrared range of frequencies, Ml 3> T, the 
contribution of the intraband absorption at low control 
voltages V is small. This implies that the last term in the 
right-hand side of Eq. (13) is much smaller than unity. 

Considering a modulator for optical radiation with 
the wavelength A = 1537 nm (Ml ~ 0.8) eV as in 
Ref. [H, setting Ml = 0.8 eV, V = V(Ml/2^/2T) 2 = 
(ed/k)(n/v w ) 2 = 337 mV, and v = 10 13 s _1 , we obtain 
that the last term in Eq. (13) is about 0.03. Thus, in the 
case of even modestly perfect GLs, this term is relatively 
small and can be omitted. This implies that in such a 
case, the modulation is primarily due to the voltage con- 
trol of the Pauli blocking but not due to the variations 
of the intraband absorption. 

In the most realistic case Ml ^> T, the ratio of the 
intensities of output (modulated by slow varying voltage) 
and input radiation, Iq and Iqo, respectively, taking into 
account that the second and fourth terms yield the same 
contribution and omitting the third and fifth terms in 
Eq. (14), can be presented as 



= exp 
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where H is the double-GL length in the directio n of ra- 
diation propagation, and hflo — hvw\/kVo/ed cx \/Vq. 
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Equation (15) describes the variation of the output radi- 
ation I caused by relatively slow variations of the applied 
voltage V m with arbitrary swing. Using Eq. (15), one can 
estimate the modulation depth too and the extinction ra- 
tio rj for the case of relatively slow modulation and when 



V m varies from V m = to V m = V™ ax 
\T max =_(ied/k)(n/v w ) 2 : mo = 1 - 



exp(- 



, where 
■p n H) and 



IV. SMALL SIGNAL LINEAR MODULATION 
CHARACTERISTICS 

We now assume that Vo is sufficiently large to form 
the degenerate electron and hole systems in the pertinent 
GLs, while the time dependence of SV m (t) <C Vq is still 
characterized by a small modulation frequency. In this 
case, from Eq. (15) we obtain the following formulas of 
the modulation amplitude dlo and the modulation depth 
of the output radiation SrriQ = SIq/I: 
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and 
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where M7o — h vw \/kVo/ed cx VM) • 

In particular, at eVo — Ml, when 5mo reaches a maxi- 
mum, Eqs. (16) and (17) yield respectively 
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V. PLASMA OSCILLATIONS IN DOUBLE-GL 
STRUCTURES 

In the double-GL structures with sufficiently high con- 
ductivity at low modulation frequencies, one can use 
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FIG. 2: Spatial distributions of the normalized amplitudes 
of the ac potential |c5<^+|/c)l/ m and \8<p-\/8V m in upper and 
lower GLs, respectively, for different modulation frequencies 
u/2tt. 



Eq. (6), i.e., put <p± — ±V m /2. However, at elevated 
modulation frequencies (for instance, in the THz range), 
the spatial distributions of the ac components of the elec- 
tron and hole charges, the electron and hole Fermi ener- 
gies, and the self-consistent electric potential are nonuni- 
form because these distributions do not follow fast modu- 
lation signals. In sufficiently perfect GLs, the modulation 
signals can excite the electron- hole plasma oscillations. 
In this situation, the ratio of the ac component of the 
radiation intensity 51^ to the input intensity is given by 
the following equation, which replaces Eq. (16): 
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(20) 

The last factor in the right-hand side of Eq. (20) accounts 
for the contributions of different parts of GLs being de- 
pendent on local values of the Fermi energies and, hence, 
the local values of the ac potentials [see Fig. 1(b)]. 

To find the distributions of 5tp+ and Stp-, one can use 
a system of hydrodynamic equations (Euler equation and 
continuity equation) adjusted to the features of the elec- 
tron and hole spectra in GLs [35j]) for the electron and 
hole plasmas in the upper and lower GLs coupled with 
the Poisson equation. For simplicity, we use the Poisson 
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equation in the gradual-channel approximation, which 
leads to Eq. (6) above. The linearized versions of the 
equations in question can be reduced to the following 
equation for the ac component of the potential at the 
frequency u> [HI]: 



d 2 Sip + lu(lu + iv) 
dx 2 s 2 



d 2 5tp- ui(u> + iv) 
dx 2 s 2 



(8<p + -S<p-) = 0, (21) 



(Sip- - 8tp+) = 0. (22) 



Here v ~ vq is the collision frequency of electrons in 
GLs with impurities and acoustic phonons and s is the 
characteristic velocity of plasma waves in GLs. Since 
electrons and holes belong to different GLs separated 
by the rather high and thick barrier, their mutual col- 
lisions can be neglected. The plasma-wave velocity is 
determined by the net dc electron and hole density (i.e., 
by the Fermi energy) E — (n/hvw) 2 /tt and the gate 
layer thickness d [3l|, HI, HE} : 



~ v w y ^aGL/iod/hvw oc 
Vq 17 ' 4 ^ 1 / 4 , where aoL = e 2 /khvw is the coupling con- 
stant (acL/& — c/vw — 300). Similar equations were 
obtained for the gated electron channels in the tradi- 
tional heterostructures. The difference is in the existence 
of two interacting channel and in different values of s . 
At d = 10 nm and hflo = 0.2 — 0.8 eV, one obtains 
s/vw ~ 3.75 — 7.5. 

The plasma wave velocity in GL system is much higher 
than that in the standard gated channels [H, HH . Large 
values of s in the double-GL should allow us to achieve 
rather high plasma frequencies (in the THz range) in the 
devices with relatively large values of L (in the microme- 
ter range) . If the left-side contact and the upper GL and 
the right-side contact and the lower GL, see Fig. 1(a)] are 
are very close [321 ] . one can disregard the gaps between 
GLs and the contacts and use the following boundary 
conditions for Eqs. (21) and (22): 
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The latter boundary condition reflects the fact that the 
electron and hole currents are equal to zero at the dis- 
connected edges of GLs (at x — — L in the upper GL and 
at x = L in the lower GL) 

Fist from Eqs. (21) and (22) we obtain 



8tp+ + 5<p- = Ax, 



(25) 



where A is a constant. Considering Eq. (25), Eqs. (21) 
and (22) can be presented as 



d 2 6<p 
dx 2 



d 2 5(p_ 
dx 2 
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(27) 



Solving Eqs. (26) and (27) with boundary conditions 
(24) and (25), we obtain 



COSJu,X 

SV m f 7„ sinj u L 



cos 7 W L 



L 



(28) 



Stp- 



SV m ( 7^ sin 7^1/ 



cos7 w L 
7^ sin 7a, L 



L 



(29) 



Here 7^ = ^/2o;(cli + iv)/ s. Introducing the charac- 
teristic plasma frequency u p = tt one obtains 
7 U L = iry/w(u) + iv)/2uj p . 

In the low modulation frequency limit (when u <C 
L 2 /2s 2 v = 4uj 2 p /tt 2 v 



-1 

' M ' 



where r M is the Maxwell 



relaxation time, and, hence, 17^1-k <C 1), one obtains 
5ip + = 8V m /2 and 8(p- — —8V m /2, i.e., the spatial dis- 
tribution of the ac potential across the GLs is flat. 

Figure 2 shows examples of the spatial distributions 
of the amplitudes of the ac potential across GLs, calcu- 
lated for different modulation frequencies using Eqs. (28) 
and (29). As seen from Fig. 2, the amplitudes |<V+| 
and \8if-\ are close to the amplitude of the applied ac 
voltage 8V m /2 at relatively low modulation frequencies 
(uj/2tt = 0.20 and 0.5 THz). However, with increasing 
u), the amplitudes dramatically increase (see the curves 
corresponding to uj/2tt = 0.80 and 1.00 THz). This is 
attributed to the excitation of plasma oscillations whose 
amplitude grows as the modulation frequency approaches 
^'^'the plasma resonance frequency (see below). 



VI. RESONANT MODULATION 

Substituting 8<p+ and 8<p- from Eqs. (28) and (29) into 
Eq. (21) and integrating over dx, we obtain 
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yielding for the normalized modulation depth: 



8m u 
5 m 



sin(7^L) 



(7^L)[cos(7 WJ L) - (7 W £) sin(7 w £)] 



(31) 



Using Eq. (31), the normalized modulation depth can 
also be expressed via the characteristic plasma fre- 
quency u p : 



2uj v 
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FIG. 3: Normalized modulation depth Sm^/Smo versus mod- 
ulation frequency ui/2n for different electron and hole collision 
frequencies v. 



FIG. 4: Normalized nodulation depth Sm^/dma versus mod- 
ulation frequency u)/2n for different plasma frequencies uj p : 
solid lines correspond to v = 1 x 10 12 s _1 , dashed lines cor- 
respond to v = 2.5 x 10 12 s _1 . 



At cu < v,u> p , Eq. (32) yields 



8m u 



tan y/: 



10JTM 



(33) 



Equations (19) and (33) provide the dependences of the 
modulation depth on the structural parameters and the 
modulation frequency (particularly the frequency depen- 
denceat the roll-off) obtained experimentally and de- 
scribed in Ref. 19. 

The character of the dependence of the modulation 
depth Sm^j on the modulation frequency oj depends on 
the collision frequency v. Figure 3 shows the frequency 
dependences (on / = u/2ir) of the normalized modula- 
tion depth Sm^/Smo calculated using Eq. (32) for devices 
with different values of the collision frequency v. It is 
assumed that the plasma frequency lu p /2it = 1.77 THz. 
For the characteristic plasma velocities s — (3.75 — 7.5) x 
10 s cm/s (as in the estimate in the previous section), this 
frequency corresponds to 2L — 0.75 — 1.5 /jm, i.e., to the 
length of GLs (size of the waveguide) close to the optical 
wavelength under consideration. 

At relatively large collision frequencies v, the modu- 



lation depth Sm^ monotonically decreases with increas- 
ing the modulation frequency [in line with Eq. (33)]. 
The curve for v = 10 X 10 12 s _1 in Fig. 3 corresponds 
to tm - 2 x 10 -13 s and the 3dB roll-off frequency 
fzdB — 0.8 THz. In the devices with a smaller lj p , i.e., 
with a smaller plasma- wave velocity, s, or larger length of 
GLs, 2L, the Maxwell relaxation time is longer, so that 
the roll-off frequency is smaller. However, in the devices 
with relatively small v <C oj p , Sm^ can exhibit a steep 
increase when oj approaches to the plasma resonance fre- 
quency. The sharpness of the modulation depth peak and 
its height rise with decreasing v (i.e., with an increasing 
quality factor of the plasma resonances, Q ~ lo p /v). This 
also seen in Fig. 4 (compare the solid and dashed lines). 

The plasma resonance frequencies, u) n oc u) p . Here 
n = 1, 2, 3, ... is the resonance index. As can be derived 
from Eq. (32), the resonance frequencies are given by the 
solution of the following equation: 



COt(7TU; n /2Wp) = (TTLU n /2u! p ). 



(34) 



Equation (34), in particular, yields iriui/2ujp ~ 0.86 < 
1, so that at uj p /2tv1.77 THz, one obtains lui/2tt ~ 
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0.97 THz. The second (even) resonance corresponds to 
w 2 /27r > 3.54 THz. Thus, apart from a pronounced first 
resonance, a fairly weak second resonance is also seen in 
Fig. 3. A characteristic plasma frequency falls in the THz 
range. As seen in Fig. 4, an increase in u p and, con- 
sequently, in cu n shifts the positions of the peaks in the 
frequency dependence of the modulation depth. How- 
ever, this shift can not be used for the electrical control 
because the bias voltage Vq corresponds to the energy 
of photons of modulated radiation. This is in contrast 
to other THz voltage-controlled devices using the reso- 
nant excitation of plasma oscillations (see, for instance, 
Ref. II). 

VII. CONCLUSIONS 

In summary, we developed a device model for an op- 
tical modulator based on the double-GL structure that 
was recently proposed and experimentally realized [l9| . 
The double-GL modulator utilizes the variation of ab- 
sorption due to the electrically controlled Pauli blocking 
effect. Our model accounts for the interband and intra- 
band absorption and the plasma effects in GLs deter- 
mining the spatio-temporal distributions of the electron 
and hole densities and the absorption coefficient, The de- 
veloped model yields the dependence of the modulation 
depth on the control voltage for strong but slow modu- 
lation signals and for small-signal modulation in a wide 



range of frequencies. The dependence of the modula- 
tion depth on the modulation frequency is determined by 
the relationship between the collision frequency of elec- 
trons and holes and the characteristic plasma frequency 
(or between the latter and the Maxwell relaxation time). 
At relatively large collision frequencies (or small plasma 
frequencies), the modulation depth is a monotonically 
decreasing function of the modulation frequency. The 
obtained dependencies qualitatively explain the experi- 
mental results. However, we predict that in the double- 
GL structures with relatively weak disorder and, hence, 
with a low collision frequencies, and a sufficiently high 
quality factor of the plasma oscillations, the modulation 
depth exhibits a sharp maximum at the modulation fre- 
quency, which corresponds to the plasma resonance. The 
frequency of the latter falls in the THz range for typical 
parameter values. This opens up the possibility to use 
the double-GL structures for effective modulation of op- 
tical radiation by THz signals. 
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